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We first review evidence for the Cooper pair insulator (CPI) phase in amorphous nanohoneycomb 
(NHC) films. We then extend our analysis of superconducting islands induced by film thickness 
variations in NHC films to examine the evolution of island sizes through the magnetic field-driven 
SIT. Finally, using the islanding picture, we present a plausible model for the appearance and 
behavior of the CPI phase in amorphous NHC films. 



INTRODUCTION 



The appearance of a giant magnetoresistance peak 
near the superconductor to insulator transition (SIT) in 
several thin film systems [iHS] has garnered a lot of recent 
attention. [H |6l [7] The factors that drive its appearance 
are not yet understood. Of the thin film systems that ex- 
hibit this exotic behavior, nanohoneycomb (NHC) amor- 
phous Bi thin films offer a unique morphology that might 
lead to an explanation. 

The amorphous NHC thin films considered here were 
fabricated by quench condensing Sb and then Bi onto 
a substrate of Anodized Aluminum Oxide (AAO) (see 
Fig. 2a). In the zero-temperature limit, an insulator- 
superconductor transition can be driven by increasing the 
NHC film thickness, or a superconductor-insulator transi- 
tion can be driven by applying a magnetic field to a super- 
conducting film. Both of these types of transitions bear 
little resemblance to their counterparts in films grown on 
planar glass substrates, even though we can be sure that 
films grown on AAO are locally homogeneous. [8 Most 
interestingly, the insulating phase in NHC films contains 
localized Cooper pairs, as confirmed by the observation of 
Little-Parks like oscillations in the magnetoresistance. [9 
Recently, the localization of Cooper pairs has been at- 
tributed to the formation of superconducting islands near 
the SIT that develop due to nano-scale undulations in 
the film thickness. [8 The appearance of these localized 
Cooper pairs is likely responsible for the puzzling features 
in the transport, including the giant magnetoresistance 
peak appearing at fields beyond the Little-Parks oscilla- 
tions (1-3T).[10^ 

In this paper we examine the transport behavior of 
amorphous NHC films near the SIT within the context 
of this Cooper pair islanding model. We present images 
indicating how increasing film thickness or magnetic field 
changes the island size. Inspired by the somewhat regu- 
lar arrangement of the islands in the images, we use an 
ordered grain array model to qualitatively explain how 
the evolution of island sizes can lead to the spectacular 
transport behavior in the Cooper pair insulator phase of 
these films. 



AMORPHOUS NHC FILM TRANSPORT NEAR 
THE SIT 

NHC thin films on the insulating side of the SIT 
show simply activated behavior, where i?n(T) = 
Roexp{To/T)^ whether they are driven by disorder or an 
applied magnetic field. Additionally, the activation en- 
ergy. To, goes continuously to zero at the transition in 
both cases. Fig la shows a typical thickness-driven in- 
sulator to superconductor transition for NHC films on 
Arrhenius plot. The linearly decreasing activation ener- 
gies of film sequences from four different experiments are 
shown in Fig lb (triangles are Tq data extracted from 
the series in Fig la). In this figure, the transition to 
superconductivity is marked by d^^^{To = 0). 

Complementary data for the field-driven 
superconductor-insulator transition are shown in 
Figs. 2c and d. The activation energy (Fig. 2d) changes 
dramatically in magnetic field, mimicking the behavior 
of the magnetoresistance with both oscillations and a 
peak in To.|9l[TO] In contrast, the prefactor, Rq^ increases 
monotonically by a factor of 2 through the entire field 
range (not shown). [10 Because of the oscillations, the 
initially superconducting film undergoes 7 consecutive 
transitions between superconducting and insulating 
states (marked by crossing the Tq = line). The R{T) 
curves transversing the SIT in c are for the field values 
chosen at = nHu to avoid the intra-oscillation 
transitions (open circles in d). The form of Tq{H) 
suggests two separate contributions to the activation 
energy: 

T^{H) = T^iH) + ToP^"^(i/). (1) 

The origin of the overarching p^^^^ feature is not yet 
determined, but a similarly shaped feature appears in 
the other highly-disordered (unpatterned) film systems, 
mentioned in the introduction and in more detail in the 
discussion. 

At low fields, we expect that Tq^^ is dominated by a 
Josephson energy between islands, which oscillates with 
the matching field period, as in Josephson Junction Ar- 
rays (JJAs).[ll The following expression is often used to 
characterize this oscillating energy scale, 

P^'^H) (x< Ejocosi^, - - 27rA,j/^o) >, (2) 
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Figure 1: a) Thickness-driven 1ST in NHC films on an Arrhenius plot. From top to bottom, the deposited film thicknesses are 
0.91, 0.92, 0.96, 0.99, 1.01, 1.03, 1.05, 1.08, 1.12 nm. b) Activation energies, obtained from the slopes of Arrhenius fits, versus 
deposited thickness for four different series of NHC films. Triangles represent To extracted from data in (a), c) Perpendicular 
field-driven SIT of a 1.13 nm NHC film on an Arrhenius plot. Field values from bottom to top are 0, 0.22, 0.45, 0.68, 0.90, 
1.11, 1.33, 1.55, 1.78, 2.00, 2.22 T. d) Activation energy versus applied field for the film in (c). The R{T) of the open circles (at 
H = uHm, Hm = 0.22 T) are shown in (c). Negative energy values are plotted for consistency; a dashed line divides negative 
and positive slopes (superconducting and insulating states). 



where Ejq is the zero-field phase coupling energy, the (/)S 
are phases on neighboring islands in the film, Aij is line 
integral of the field's vector potential between islands, 
and ^0 = h/{2e) is the magnetic flux quantum. The 
presence of this term in the activation energy is strong 



evidence that activated Cooper pair tunneling dominates 
the transport. 
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COOPER PAIR ISLANDING IN NHC FILMS 

A comparison of NHC and uniform films with increas- 
ing deposition thickness intimates that structure in the 
AAO substrates leads to Cooper pair localization. The 
A AO is not totally flat: regular height variations appear 
around the holes, as can be seen in the atomic force mi- 
croscope image in Fig. 2a. Since evaporated material 
impinges normal to the (average) plane of the substrate, 
thickness variations develop according to the local slope: 

d(x, y) = d^^P ^ = , (3) 

where d{x^ y) is the local film thickness and h{x^ y) is the 
local height of the substrate. Evidently, the film grows 
thickest in the flat regions of the substrate and thinnest 
in the steep regions. Fig. 2b shows the height profile 
along a line scan of the substrate as well as the corre- 
sponding thickness profile, calculated using Eqn. 1 and 
^dep _ Q ggj^m (an insulating film). Also on this plot, 
we mark the critical thicknesses for conduction, d^l^ and 
superconduction, (ifgx, in reference films. It is evident 
that some regions of the film are thick enough to support 
superconductivity while others are not. 

By applying Eqn. 1 to the AFM image of Fig 2a, we 
create a map of the local film thickness for a given d^^P. 
In addition, we color the map so that regions of film 
that conduct (c^q^ < d{x,y) < <ifsT) pink, and re- 
gions that superconduct {d{x^y) > <ijgx) are blue. Non- 
conducting film {d{x^y) < dQ^.) is gray. Fig. 2c shows 
the resulting morphology of an insulating film near the 
1ST {d^^^ = 0.89nm). Here, regions of superconductivity 
are interspersed in regions of simply conducting film. 

Taking the analysis one step further, we can transform 
the film thickness variation map into a coupling constant 
map, again by comparing to data of reference films. In 
uniform superconducting films, the transition tempera- 
ture depends on film thickness as. 

For a-Bi, T^^^^^ = 6 K and dgrj. = 0.7nm.(Jay thesis) 
Assuming this relation to be true on a microscopic level, 
we add a color bar to the right side of Fig. 2b to con- 
vert our map of local film thicknesses to one of local TcS. 
Considering the picture as a whole, the film appears to 
be made up of an array of superconducting islands con- 
nected by non-superconducting film. The local coupling 
constant, marked by a local Tc, varies over an island. It 
is strongest in the center, where d{x,y) is largest, and 
weakest on the edges. By analyzing these maps through 
the 1ST, we find that the superconducting islands grow 
with (i^^P and coalesce at the transition. [8 

Viewing the transitions as the coalescence of conduct- 
ing or superconducting regions of film allows us to under- 
stand the critical deposition thicknesses for conduction 



{d^f^) and superconduction {d^§^) in NHC films.[8l 
This result is corroborated by modeling the film as a 
network of islands linked by resistor elements, shown in 
Fig 2d. At the SIT, the resistance of a single weak link 
in this network is consistently close to Rq = /i/4e^, the 
critical value of sheet resistance in uniform films. Both of 
these observations indicate that the weak links between 
superconducting islands control the transition to super- 
conductivity in these films. [8] 

SHRINKING COOPER PAIR ISLANDS WITH 
AN APPLIED MAGNETIC FIELD 

We can extend our islanding picture to undulating 
films in an applied magnetic field by once again compar- 
ing the local film thicknesses to uniform films. This fairly 
primitive approach has led to a nice picture that quali- 
tatively agrees with some models described in the dis- 
cussion. We attribute the critical field in uniform amor- 
phous Bi films of various thicknesses, Hc2{d^^^), to local 
film thicknesses in undulating amorphous Bi films. In 
uniform films, Hc2 depends on thickness as, 

Hc2 = 7.12T(1 - 0.7nm/^^). (5) 

(determined by a fit to the data)[T2^. Using this phase 
line, we can define a critical thickness for superconduc- 
tivity for every value of applied field, dsiriH). In Fig. 4, 
the superconducting film of Fig. lc,d (d^^P = 1.08nm) is 
shown for a series of applied magnetic fields. At each field 
value the color bar is redefined to include a region of pre- 
viously superconducting film made normal by the applied 
field (cranberry colored). In the figure, the islands shrink 
with increasing field, eventually disappearing at H=2.5T 
(not shown). The accompanying histograms show the 
magnetic field evolution of the island size distribution, as 
the number of islands of each size (in a Ixl/im region), 
and the percent film area taken up by islands of that 
size. At H=OT, a single island dominates the film area, 
but as the field is increased these islands separate. The 
second panel from the top shows the rather broad dis- 
tribution of islands at the field-driven SIT (0.8T), where 
they have just separated. At i^peak=2T, the islands are 
much smaller and fairly normally distributed in number 
and area around the estimated coherence area, tt^^. Be- 
yond the peak, at H=2.4T, the islands are all smaller 
than TT^^ and have nearly disappeared. 

DISCUSSION 

Dramatic transport behavior near the SIT, similar to 
what has been presented here (see Fig. 1), has been 
observed in a number of thin film systems, including in- 
dium oxide, titanium nitride, and beryllium, [l] |4l US] In 
particular, a number of models have been formulated 
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Figure 2: Figure adapted from reference [8^ a) AFM image of the A AO substrate. The (white) scale bar spans 100 nm. 
b) Height profile of substrate along a line scan (black line in (a)). The dashed gray curve is the corresponding film thickness 
profile calculated using Eqn. 3. Pink and blue dashed lines mark the critical thicknesses for conductivity and superconductivity, 
respectively, c) Film thickness variations calculated using the image in (a) and Eqn. 3. Colored according to critical thicknesses 
in a reference film: pink if the film is thick enough to conduct (d^^^ > (TqI = 0.5nm), and blue if the film is thick enough to 
superconduct (d^^^ > c/iIt = 0.7nm). The right color bar defines the local coupling constant, Tc, calculated using Eqn. 4. d) 
A model of the film as superconducting islands connected by resistor elements, or weak links. 



to try to explain the origin of the giant magnetoresis- 
tance peak. [3l |6l [3 [14] In these models, the peak is com- 
monly attributed to the existence of islands of incoherent 
Cooper pairs. Our analysis of the morphology of NHC 
films suggests a well defined mechanism for island forma- 
tion. This mechanism provides a picture of the geomet- 
rical arrangement and rough sizes of the islands through 
the disorder and magnetic field tuned SITs. Below, we 
combine these unique insights with and ordered grain ar- 
ray model to provide a qualitative explanation for the 
activated insulator and MR peak. 

The hard gap at the SIT 

A number of models produce energy scales related to 
Cooper pairs in the insulator phase that may correspond 
to the activation energy seen in insulating NHC films (see 
Fig. lb,d). Trivedi and coworkers[15^ showed that an 
attractive Hubbard model with disorder predicts a pair 
gap that opens at the SIT and grows into the insulating 
state. This gap corresponds to the energy required to 



inject a Cooper pair into the insulator. Very recently, 
Miillerp^ [TB] proposed that a mobility edge can appear 
in a disordered, localized bosonic system with Coulomb 
interactions. The distance to the mobility edge, which 
serves as the activation energy for Cooper pair trans- 
port, increases with disorder and magnetic field. Finally, 
and of primary interest here, an ordered array of small, 
Josephson coupled, superconducting grains is expected 
to exhibit a Mott insulator phase with an energy barrier 
for Cooper pair injection (see [17j and refs. therein). In 
the limit of very weak coupling between grains this bar- 
rier corresponds to the self charging energy for a grain. 
Close to the SIT, the increased coupling between 
grains, as measured by the Josephson energy, Ej^ leads 
to a downward renormalization of this barrier, Eq Eq. 
Deep in the insulator, Eq :^ Ec — E:^/A. As the intergrain 
normal state resistance approaches i?Q, this renormaliza- 
tion can be very substantial, (citation) 

The network of islands that results from CPs local- 
ized by coupling constant variations suggests applying 
the ordered grain array model to NHC films. Within 
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Figure 3: Superconducting islands in an increasing magnetic field for the l.OSnm film of Fig lc,d. Right: Spatial map of SC 
islands (blue) in an insulating background (pink) with increasing magnetic field. The cranberry color represents regions of 
initially superconducting film for which Hc2{d{x,y)) < i^appiied, calculated using Eqn. 5. Left: Distribution of superconducting 
island size and number of islands of each size for the applied field series. The first, third, and fourth panels (from the top) 
correspond to R(T) shown in Fig. Ic. 



this model, we estimate that the islands contain few 
Cooper pairs and have a large energy level spacing. The 
bare charging energy between grains, Eq^ is found to be 
enormous for the islands in Fig. 2 (approximately 100 
meV)[9 and much greater than the Josephson energy, 
Ej^ so that the applicability of the model would seem to 
rely on a very strong renormalization of Eq near the SIT. 

The number of Cooper pairs in a coherence volume can 
be estimated as A^cp = island ^(^f) A, where Fisiand is the 
volume of the superconducting island, t^iep) is the den- 
sity of states at the Fermi energy and A is the amplitude 
of the superconducting order parameter. Taking an is- 



land near the SIT to be 20nm in diameter and Inm thick, 
^island ^ lOOTTum^. We estimate iy{eF) for a-Bi to be 
2 X 10^^cm~^eV~^, the value for crystalline Pb. Finally, 
we estimate A using 2 A = S.bksTc and the local maxi- 
mum Tc from the coupling constant map (Fig. 2c). Us- 
ing these numbers we find that Ncp ^ 2. Additionally, 
the energy level spacing in an island, S = l/i^(ei?)l^siand7 
yields S ^ 1.76K. 

To estimate the Josephson energy, we use Ej = 7r^A/2 
where the normalized conductance, ^ ~ 1 near the tran- 
sition. For Tc = 2K, this expression gives Ej ^ 5K(^ 
0.4meV). Finally, we calculate the charging energy with 
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Ec = 4e^/eoeL using L = 20nm, as the superconducting 
island size (see Fig. 2c). We find Eq ^ lO^K/e, implying 
Ec ^ ^, ^c, Ej for any reasonable e. 

Near the SIT, the condensation energy per grain, the 
Josephson energy, the energy level spacing, and Tc are 
all of the same order. The charging energy per grain 
is very large, however. We expect it to be dominant in 
this model even after the necessary renormalization. On 
approaching the disorder-driven SIT from the insulating 
side, the superconducting islands grow with d^^^ and co- 
alesce at the transition. [8 As the islands grow, their sep- 
aration shrinks and they become better coupled. These 
changes result in a diminishing Eq and an increasing Ej 
as d^^^ (^^Ft- The growing islands accommodate in- 
creasingly more Cooper pairs. These changes result in 
a reduction of Eq, which corresponds to the decreasing 
To(^i^^P) in Fig. lb. 



The grain array model attributes the activation energy 
to charging effects, which other models do not need to in- 
clude (see [H]). It also requires the array to be ordered, 
as disorder in the array would likely lead to variable range 
hopping or percolation phenomena. The model of Dubi 
and Meir[6] invokes magnetic field induced shrinkage of 
superconducting islands, but also uses a percolation con- 
duction model because many of the systems exhibiting 
the MR peak are disordered arrays. 

We hope our results provide a starting point for ex- 
plaining the common features in the transport across 
seemingly dissimilar systems exhibiting Cooper pair insu- 
lating behavior. Future theoretical work needs to address 
the role played by disorder and the absence of quasiparti- 
cle transport in these Cooper pair insulators. Further ex- 
periments that can establish the importance of Coulomb 
interactions and disorder to these SITs are also necessary. 



The magnetoresistance peak 

Here, we extend the discussion of NHC films modeled 
as an ordered island array to account for the effects of an 
applied magnetic field on their transport. With Eqn. 2, 
we attribute the influence of the magnetic field on the ac- 
tivation energy to a combination of two terms. The low 
field magneto-oscillations in Tq®^ (Fig. Id) result from 
field-induced modulations of Ej and, as a result, Eq- At 
higher fields, Tq^^^ dominates and gives rise to the mag- 
netoresistance peak. j'P®^^ grows due to the magnetic 
field induced shrinking of the islands, which drives up Ec 
and reduces the inter-island coupling, reducing Ej. Also, 
the magnetic field reduces A, which also reduces Ej. Al- 
together, these effects drive up Eq, or Tq^^^. Eventu- 
ally, the magnetic field depresses A below Eq, at which 
point quasiparticle tunneling rather than CP tunneling 
becomes energetically favorable. As pointed out by Gant- 
makher [3 and others, the magnetoresistance becomes 
negative for quasiparticle tunneling because the tunnel- 
ing rate increases as the field depresses the gap. 

Conclusion 

We reviewed evidence for the existence of a Cooper pair 
insulator phase in amorphous NHC films. AFM images 
reveal that structure in the substrates induce the forma- 
tion of superconducting islands separated by weak links. 
We presented the evolution of these island sizes with film 
thickness and magnetic field. We then discussed an or- 
dered grain array model that can qualitatively account 
for the thickness and magnetic field tuned SITs as well as 
the large magnetoresistance peak observed in the Cooper 
pair insulator state. 

While the ordered island array scenario is plausible, ex- 
perimental observations still allow other interpretations. 
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